CLEO CONF 00-1 
ICHEP00-78 



B — > D*7T + 7r 7T 7T°, D^cjiT and the Observation of a 
Wide 1~ cj7r~ Enhancement at 1418 MeV 

CLEO Collaboration 
June 15, 2000 



Abstract 

We report on the observation of B — > D* n + n~ n~ ir° decays. The branching ra- 
tios for D* + and D*° are (1.72±0.14±0.24)% and (1.80±0.24±0.25)%, respectively. 
Each final state has a D*ujit~ component, with branching ratios (0.29±0.03±0.04)% 
and (0.45±0.10±0.07)% for the D* + and D*° modes, respectively. We also observe 
B — > Dljtt~ decays. The branching ratios for D + and D° are (0.28±0.05±0.03)% and 
(0.41±0.07±0.04)%, respectively. The uji:~ appears to come from the decay of a wide 

resonance. A fit to a Breit-Wigner shape gives a mass of 1418±26±19 MeV and 
width of 388±44±32 MeV. We identify this object as the p(1450) or pf . 
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1 Introduction 



Understanding hadronic decays of the B is crucial to insuring that decay modes used for measure- 
ment of CP violation truly reflect the underlying quark decay mechanisms expected theoretically. 

Currently, measured exclusive branching ratios for hadronic B decays total only a small fraction 
of the hadronic width. The semileptonic branching ratio for B — * Xe~v, Xfi~v, and Xt~v totals 
approximately 25% Q. The measured hadronic decay modes for the B including D + (mr~), 
D* + (nir~), where 3 > n > 1, D + ^Ds { *\ and J/i/j exclusive totals only about 10% Q. (The B~ 
modes total about 12%.) Thus our understanding of hadronic B decay modes is not yet well based 
in data. 

It is also interesting to note that the average charged multiplicity in a hadronic B° decay is 
5.3±0.1 Since this multiplicity contains contributions from the D + or D* + normally present 
in B decay, we expect a sizeable, approximately several percent, decay rate into final states with 
four pions ||. The seen D^*\mr)~ final states for n < 3 are consistent with being quasi-two-body 
final states. For n of two the p~ dominates, while for n of three the aj~ dominates Qj. These decays 
appear to occur from a simple spectator mechanism where the virtual W~ materializes as a single 
hadron: 7r~, p~ or . 

In this paper we investigate final states for n of 4. We will show a large signal for the 
D* + ir + ir~ ir~ ir° final state in section ||. In section || we will show that a substantial fraction, 
~20% arise from D* + uir~ decays and that the ujir~ mass distribution has a resonant structure 
around 1.42 GeV with a width of about 0.4 GeV. In section [B|, the similar conclusions are drawn 
about the .D*°7r + 7r~7r~7r° final state. The same structure is shown to exist in Dojit~ final states 
(section |?|) and we will use these events to show in section ffl] that the spin-parity is most likely 1~. 
This state is identified as the p', also sometimes called the /o(1450). Other resonant substructure is 
searched for, but not found (section ^). Finally we summarize our findings and compare with the 
predictions of factorization and other models in section [!(]. 

The data sample consists of 9.0 fb _1 of integrated luminosity taken with the CLEO II and II. V 
detectors Q using the CESR e + e~ storage ring on the peak of the T(4S) resonance and 4.4 fb _1 in 
the continuum at 60 MeV less center-of-mass energy. The sample contains 19.4 million B mesons. 

2 Common Selection Criteria 

Hadronic events are selected by a minimum of five charged tracks, total visible energy greater 
than 15% of the center-of-mass energy, and a charged track vertex consistent with the nominal 
interaction point. To reject continuum we require that the Fox- Wolfram moment R2 be less than 
0-3 §. 

Track candidates are required to pass through a common spatial point defined by origin of all 
tracks. Tracks with momentum below 900 MeV/c are required to have ionization loss in the drift 
chamber within 3<7 of their assigned massj^j (These requirements are not imposed on slow charged 
pions from D* + decay.) Photon candidates are required to be within the "good barrel region," 
within 45° of the normal to the beam line, and have an energy distribution in the Csl calorimeter 
consistent with that of an electromagnetic shower. To select 7r°'s, we require that the diphoton 
invariant mass be between -3.0 to +2. 5a, where a varies with momentum and has an average value 
of approximately 5.5 MeV. For each two-photon mass combination a is calculated. After candidate 
selection the two-photon's are kinematically fit by constraining their invariant mass to that of the 

1 Here and throughout this paper a indicates an r.m.s. error. 
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We select D° and D + candidates via the decay modes shown in Table |l|. We require that the 
invariant mass of the D candidates lie within ±2. 5a of the known D masses. The ex's are also listed 
in Table |l|. The D° widths vary with the D° momentum, p, (units of MeV). 

We select D* + candidates by imposing the addition requirement that the mass difference be- 
tween tt + D° and D° combinations is within ±2.5<r of the known mass difference. For the D*°, we 
use the same requirement for the ir°D° decay. The mass difference resolutions are 0.63 and 0.90 
MeV, for the ir + D° and 7r°D° modes, respectively^] 



Table 1: Mass Resolutions (a) in MeV 



D+ -> K-Ji + Ti + 


D° -> K~7r + 


D° - 


-> K~Tl + TT° 


D° -> K-n+n+n- 


6.0 


px0.93xl0- 3 +6.0 


px0.6i 


}xl0- 3 +11.6 


px0.92xl0~ 3 +4.7 



3 Observation of B — > D* + 7r + 7r tt tt° Decays 

3.1 B Candidate Selection 

We start by investigating the D* + (47r)~ final state.^ The D* + candidates are combined with all 
combinations of 7r + 7r~7r~7r° mesons. 

Next, we calculate the difference between the beam energy, Ei, earn , and the measured energy 
of the five particles, AE. The "beam constrained" invariant mass of the B candidates, Mg, is 
computed from the formula 

M 2 B = E 2 beam -(Y,p!) 2 ■ (1) 

i 

To further reduce backgrounds we define 

2 / AM D * \ 2 f AM D \ 2 v / AM n o \ 2 
Xb U(AM D ,)/ U(AM D )j + ^ U(AM l0 )j ' {> 

where AMd* is the computed D* — D° mass difference minus the nominal value, AMp is the 
invariant candidate D° mass minus the known D° mass and AM n o is the measured 77 invariant 
mass minus the known tt° mass. All 7r°'s in the final state are included in the sum. The cr's are 
the measurement errors. We select candidate events in each mode requiring that xt < Cn, where 
C n varies for each decay D° decay mode. For the Knit decay modes we use C n = 12, 8, and 6, 
respectively. 

3.2 Branching Fraction and (4tt) _ Mass Spectrum 

We start with the D° — > K~ir + decay mode. We show the candidate B mass distribution, Mb, 
for A.E in the side-bands from -5.0 to -3. Oct and 5.0 to 3. 0<r on Fig. [l](a). The AE resolution is 
18 MeV (a). This gives a good representation of the background in the signal region. We fit this 
distribution with a shape given as 

back(r) = Pl rVl^e- p2{1 - r2) , (3) 

where r = Mb/5.2895, and the pi are parameters given by the fit. 

2 In this paper (4-7r)~ will always denote the specific combination 7r + 7r~7r~7r°. 
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Figure 1: The B candidate mass spectra for the final state Z}* + 7r + 7r~7r~7r°, with D° — > K~tt + 
(a) for AE sidebands and (b) for AE consistent with zero. The curve in (a) is a fit to the 
background distribution described in the text, while in (b) the shape from (a) is used with 
the normalization allowed to float and a signal Gaussian of width 2.7 MeV is added. 

We next view the Mg distribution for events having AE within 2a around zero in Fig. |l](b). 
This distribution is fit with a Gaussian signal function of width 2.7 MeV and the background 
function found above whose normalization is allowed to vary. We find 358±29 events in the signal 
peak. 

We repeat this procedure for the other two D° decay modes. The Mb spectrum for AE 
sidebands and signal region is shown in Fig. ^. The AE resolution is 22 MeV in the K~ir + ir° mode 
and 18 MeV in the K~n + n + ii~ mode. Signal to background ratios are worse in these two modes, 
but the significance is still large. The number of signal events in each mode are shown in Table ^. 

Table 2: Event numbers for the D* + 7t + 7t~7t~tt final state 



D° 


Decay Mode 


Fitted # of events (%) 


K~ 


7T + 


358±29 


K~ 




543±49 


K~ 




329±41 




We choose to determine the branching fraction using only the D° — > K~ir + decay mode because 
of the relatively large backgrounds in the other modes and the decreased systematic error due to 
having fewer particles in the final state. In order to find the branching ratio we use the Monte 
Carlo generated efficiency, shown in Fig. ||| as a function of (47r)~ mass. The efficiency falls off at 
larger (47r) _ masses because the detection of the slow ir + from the D* + decay becomes increasingly 
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Figure 2: The B candidate mass spectra for the final state -D* + 7r + 7r~7r~7r°, the left-side plots 
are for D° — > K^-k + -k° (a) AE sidebands, (b) for AE consistent with zero; the right-side 
plots are for D° — > f^~7r + 7r + 7r~ (c) A£? sidebands, (d) for AE consistent with zero. The 
curves in the top plots, (a) and (c), are fits to the background distribution described in the 
text, while in the bottom plots, (b) and (d), the shapes from (a) and (c) are used with the 
normalization allowed to float and a signal Gaussian of width 2.7 MeV is added. 

difficult. Since the efficiency varies with mass we need to determine the (4-7r) - mass spectrum. To 
rid ourselves of the problem of the background shape, we fit the B candidate mass spectrum in 
50 MeV bins of (47r)~ mass. (The mass resolution is approximately 12 MeV.) The resulting (4-7r) - 
mass spectrum is shown in Fig. |]. There are indications of a low mass structure around 1.4 GeV. 
This will be investigated further in this paper. 
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Figure 3: The efficiency for the final state D* + n + n n tt°, with D° — > K 
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Figure 4: The invariant mass spectra of n + n n n° for the final state D* + n + n n n°, with 
D° — > i^ _ 7r + , found by fitting the 5 yield in bins of Air mass. 



We find 

B(B° ^ D*+tt + it + it-it ) = (1.72 ±0.U± 0.24)% . (4) 

The systematic error arises mainly from our lack of knowledge about the tracking and tt° 
efficiencies. We assign errors of ±2.2% on the efficiency of each charged track, ±5% for the slow 
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pion from the D* + , and ±5.4% for the ir°. The error due to the background shape is evaluated in 
three ways. First of all, we change the background shape by varying the fitted parameters by la. 
This results in a change of ±3%. Secondly, we allow the shape, P2, to vary (the normalization, p±, 
was already allowed to vary). This results in 3.8% increase in the number of events. Finally, we 
choose a different background function 

back'(r) = p\r \f~V — r 2 (\+ p2r + p^r 2 + p^r 3 ^ , (5) 

and repeat the fitting procedure. This results in a 3.7% decrease in the number of events. Taking 
a conservative estimate of the systematic error due to the background shape we arrive at ±3.8%. 
We use the current particle data group values for the relevant D* + and D° branching ratios of 
(68.3±1.4)% (D* + -> ir + D°) and (3.85±0.09)% (D° -> K~it + ), respectively 0. The relative 
errors, 2.0% for the D* + branching ratio and 2.3% for the D° are added in quadrature to the other 
sources of systematic error, yielding a total systematic error of 10%. 

We wish to search for narrow structures. However, we cannot fit the B mass spectrum in small 
(47r)~ mass intervals due to a lack of statistics. Thus we plot the (47r)~ mass for events in the 
Mb peak for the D° — > K~ir + mode and the sum of all three modes in Fig. ||. We also plot 
two background samples: events at lower Mb (5.203 - 5.257 GeV) and those in the AE sideband 
separately. First we view the plots in the canonical 50 MeV bins. Both background distributions 
give a consistent if somewhat different estimates of the background shape. (Each background 
distribution has been normalized to the absolute number of background events as determined by 
the fit to the Mb distribution.) In any case no prominent narrow structures appear in the histograms 
for the 10 MeV binning. 

4 The B° -> D* + uj7r~ Reaction 

To investigate the composition of the (47r)~ final state, we now investigate the 7r + 7r _ 7r° mass 
spectrum for the events in the B peak. All three D° decay modes are used. We show the tt + tt^tt 
invariant mass distribution for events in the B mass peak in Fig. |6| (there are two combinations per 
event). A clear signal is visible at the uj. The histograms on the figure are for events either in the 
lower Mb range, from 5.203 GeV to 5.257 GeV, or in the previously defined AE sidebands; no uj 
signal is visible. 



9 



> 



O 

in 
w 

£ 
(D 
> 
LU 



40 



20 



Signal 

M B sideband 
AE sideband 




3.5 



(GeV) 



200 



150 



o 

c 

> 

LU 



100 



50 



Signal 



M B sideband 



AE sideband 




3.5 



(GeV) 




M 




7i ji Jul (GeV) WL,,, (GeV) 



Figure 5: The invariant mass spectra of 7r + 7r~7r~7r° for the final state D* + 7t + 7t~tt~'k , with 
D° — > ^^71+ (upper left), and the sum of all three D° decay modes (upper right). Events are 
selected by being within 2a of the B mass. The solid histogram is the background estimate 
from the M B lower sideband and the dashed histogram is from the AE sidebands; both are 
normalized to the fitted number of background events. The same distributions in smaller 
bins (lower plots). 
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Figure 6: The invariant mass spectra of 7r + 7r~7r° for the final state D* + n + n~ n~ n° for all 
three D° decay modes. The solid histogram is the background estimate from the Mb lower 
sideband and the dashed histogram is from the AE sidebands; both are normalized to the 
fitted number of background events. 



The purity of the u sample can be further improved by restricting candidates to certain regions 
of the Dalitz plot of the decay products. We define a cut on the Dalitz plot as follows. Let To, T + 
and T_ be the kinetic energies of the pions, and Q be the difference between the u mass, M w , and 
the mass of the 3 pions. We define two orthogonal coordinates X and Y, where 



X = 3T /Q-1 

Y = y/3(T+-T-)/Q . 

The kinematic limit that defines the Dalitz plot boundary is defined as 

1 



Y 2 

boundary 



"\^boundary ~f~ 1) {-^boundary ~t~ 1 ~f~ 0") (1 ~f~ bj (^boundary ~f~ 1 c)) 



(6) 
(7) 



(8) 



where a = Qm /Q, b = 6m 2 /(M OJ Q), c = 3(M a) - m ) 2 /(2M a; Q), m is the mass of a charged pion 
and mo the mass of the neutral pion. 

For any set of three pion kinetic energies, we define a variable r, properly scaled to the kinematic 
limit as 



\ 



X 2 + Y 2 



X 2 4- Y 2 

boundary ' boundary 



0) 



where the boundary values are found by following the radial vector from (0,0) through (X, Y). 

For events in the B mass peak we show in Fig. ^ the 7r + 7r~7r° invariant mass for three different 
cuts on r. The to signal is purified by using a selection on r. 
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Figure 7: The invariant mass spectra of 7r + 7r~7r° for the final state D* + 7t + 7t~7t~tt for all 
three D° decay modes for three selections on r: (a) 1, (b) 0.7 and (c) 0.5. The solid histogram 
is the background estimate from the Mb lower sideband and the dashed histogram is from 
the AE sidebands; both are normalized to the fitted number of background events. 



For further analysis we select oj candidates within the the tt + tt^tt mass window of 782±20 
MeV with r < 0.7. We also abandon the \ 2 cu t as background is less of a problem. In Fig. || 
we show the B candidate mass distribution for the D* + ojtt~ final state summing over all three D° 
decay modes. (The signal is fit with the same prescription as before.) There are 136±15 events in 
the peak. 

In Fig. ^ we show the ojh~ mass spectrum in the left-side plot. The solid histogram shows 
events from the lower Mb sideband region suitably normalized. The dotted histogram shows the 
background estimate from the AE sidebands, again normalized. In the signal distribution there 
is a wide structure around 1.4 GeV, that is inconsistent with background. We re-determine the 
u)ir~ mass distribution by fitting the Mb distribution in bins of luit~ mass, and this is shown on 
the right-side. Fitting to a Breit-Wigner function, we find a peak value of 1416±37 MeV and a 
width of 402±47 MeV. These numbers change to 1432±37 MeV and 376±47 MeV, respectively, 
after applying a correction for the variation of efficiency with mass. 
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Figure 8: The Mb spectra for D* + ujtt~ for all three D° decay modes, (a) AE sidebands 
and (b) AE around zero. 
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Figure 9: The invariant mass spectra of ujtt~ for the final state D* + 7t + 7t~7t~tt for all three 
D° decay modes, (left) The solid histogram is the background estimate from the Mb lower 
sideband and the dashed histogram is from the AE sidebands; both are normalized to the 
fitted number of background events, (right) The mass spectrum determined from fitting the 
M B distribution and fit to a Breit-Wigner function. 
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Knowing the ujtt mass dependence of the efficiency we evaluate the branching fraction: 

B(B° -» D* + lott-) = (0.29 ±0.03 ±0.04)% . (10) 



We tentatively label the state at 1432 MeV the A" and investigate its properties later. The 
ujir~ comprises about 17% of the (47r) _ final state. All of the u}ir~ final state is consistent with 
coming from A~ decay. 



5 Observation of B 



D*°7T + 7T-7T-7T° 



We proceed in the same manner as for the B reaction with the exception that we use the 
D*° — > tt°D° decay mode and restrict ourselves to the D° — ► K~ir + decay mode only due to 
large backgrounds in the other modes. The x 2 is calculated according to equation § and we use a 



cut value of 8. The Mb distributions for AE sidebands and signal data are shown in Fig. 10 for 
the D° — > K~ir + decay mode. We see a signal of 195±26 events yielding a branching fraction of 



B(B- 



JJ TT TT IT IT 



(1.80 ±0.24 ±0.25)% 



(11) 
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Figure 10: The B candidate mass spectra for the final state .D*°7r + 7r~7r~7r°, with D° — > 
K~tt + (a) for AE sidebands and (b) for AE consistent with zero. The curve in (a) is a fit 
to the background distribution described in the text, while in (b) the shape from (a) is used 
with the normalization allowed to float and a signal Gaussian of width 2.7 MeV is added. 



The tt + it~it mass spectrum shown in Fig. 11 shows the presence of an u>. Selecting on the 
presence of an lu with r < 0.7 we show the sideband and signal plots in Fig. (Here we do not 
use the previously defined \ 2 cut.) The branching ratio is 



D* 



UJ1T 



(0.45 ±0.10 ±0.07)% 



(12) 
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The invariant mass spectra of 7r + 7r~7r° for the final state D*°tt + ~7t~7t~7t° for the 
r + decay mode. The solid histogram is the background estimate from the Mb 
lower sideband and the dashed histogram is from the AE sidebands; both are normalized to 
the fitted number of background events. There are two combinations per event. 
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Figure 12: The M B spectra for D*°ujtt for the D° — > K n + decay mode, (a) AE sidebands 
and (b) AE around zero. 
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In Fig. ^ we show the u>ir~ mass spectrum. We see an enhancement at around 1.4 GeV as in 
the neutral B case. A fit to the data gives a mass of 1367±75 MeV and width of 439±135 MeV, 
consistent within the large errors with the B case. (We do not have enough statistics here to fit 
the Mb distribution in bins of lott~ mass.) The ujit~" fraction of the (4-7r) - final state is 25%, and 
all the ojtt~ is consistent with coming from the A~ . 
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Figure 13: The invariant mass spectra of um~ for the final state Z}*°7r + 7r~7r~7r° for the 
D° — > K~7i + decay mode, (left) The solid histogram is the background estimate from the 
Mb lower sideband and the dashed histogram is from the AE sidebands; both are normalized 
to the fitted number of background events, (right) The data fit to a Breit-Wigner signal and 
a smooth background function. 



6 Analysis of D* + ujtt Decay Angular Distributions 

The A~ is produced along with a spin-1 D* from a spin-0 B. If the A~~ is spin-0 the D* would 
be fully polarized in the (J, J z ) = (1,0) state. If the A~ were to be spin-1 any combinations of 
z-components would be allowed. It is natural then to examine the helicity angle of the D* + by 
viewing the cosine of the helicity angle of the n + with respect to the B in the D* + rest frame. 

Another decay angle that can be examined is that of the ujtt system. If the A~ is spin-0, the 
<jj is polarized in the (1,0) state and may be if the A~ is spin-1. Here the helicity angle is defined 
as the angle between the normal to the u decay plane and the direction of the A~ in the u rest 
frame. For a spin-0 A~ the distribution will be cosine-squared. Again full polarization is possible 
if the A" is other than spin-0, but any distribution other than cosine-squared would demonstrate 
that the spin is not equal to zero. 

For this analysis we use all three D° final states for the D* + final state. To find the distributions 
we fit the number of events in the Mb candidate plot selected on different angle bins. The uir mass 
is required to be between 1.1 and 1.9 GeV. This leaves 111±13 events. 
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In Fig. 14 we show the helicity angle distribution, cos for the D* decay. The data have been 



corrected for acceptance. We also show the expectation for spin-0 from the Monte Carlo. The data 
have been fit for the fraction of longitudinal polarization. We find 



r 



= 0.63 ±0.09 



(13) 




Figure 14: The cosine of the angle between the D° and the D* flight direction in the D* 
rest frame for the D*A~ final state (solid points) along with a fit (solid curve) allowing the 
amount of longitudinal and traverse polarization to vary. The dotted curve is the expectation 
for a spin-0 A". 

The cosOd* distribution is not consistent with full polarization, yielding a x 2 °f 17.7 for 5 
degrees of freedom. The helicity angle distribution for the A~ —* luit~ , cos 6^, is shown on Fig. 15. 
Furthermore the cos^ distribution is quite inconsistent with a cos 2 9 W , yielding a x 2 of 109 for 5 
degrees of freedom. Therefore, we rule out a spin-0 assignment for the A~ . 

To determine the J p we need a more well defined final state. This is provided by analysis of 
B — ► Dum~ decays. 



7 Observation of B —> Dujtt Decays 



7.1 B candidate selection 



K-TT+ 



or 



D+ 



K 7T + 7r + decay. 



Here we study the reactions B —>■ Dujit~ , with either a D° 
Other D° or D + decays have substantially larger backgrounds. 

Although we are restricting our search to u's, we define two Tr + ir~ir° samples. One within 20 
MeV of the known uj mass (782 MeV) and the other in either low mass or high mass sideband 
defined as three tt mass either between 732 and 752 MeV or between 812 and 832 MeV. We also 
require a cut on the uj Dalitz plot of r < 0.7. 
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Figure 15: The cosine of the angle between the normal to the ir + 7r~ir decay plane and 
the uj boost direction for the D*A~ final state (solid points) along with a fit (solid curve) 
allowing the amount of longitudinal and traverse polarization to vary. The dotted curve is 
the expectation for a spin-0 A~ . 



To reduce backgrounds we define 

2 = / AM D \ 2 / AM, \ 2 / AM w o 
Xb \a(AM D )) + {a(AM UJ )J \a(AM w o) 

where AM/j is the invariant candidate D° mass minus the known D° mass, AM U is the invariant 
candidate uj mass minus the known uj mass, and AM n o is the measured 77 invariant mass minus 
the known ir° mass. The cr's are the measurement errors. We select candidate events requiring that 
xl is < 12 for the Ktt mode and <6 for the Ktttt. 




7.2 B -> D°ojh- Signal 

We start with the D° — > K~ir + decay mode, for events in the uj peak. We show the candidate B 
mass distribution, Mb, for AE in the side-bands from -7.0 to -3.0a and 7.0 to 3. 0<r on Fig. |l6|(a). 
The AE resolution is 18 MeV (a). This gives a good representation of the background in the signal 
region. We fit this distribution with a shape given as 

back(r) = Pl rVl^e- p2{1 - r2) , (15) 

where r = Mb/5.2895, and the pi are parameters given by the fit. 

We next view the Mb distribution for events having AE within la around zero in Fig. |l6|(b). 
This distribution is fit with a Gaussian signal function of width 2.7 MeV and the background 
function found above whose normalization is allowed to vary. We find 88±14 events in the signal 
peak. 
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Figure 16: The B candidate mass spectra for the final state D°ujtt~ , with D° 



K-ti + . 



(a) 



for AE sidebands, and (b) for AE consistent with zero. The vertical scale in (a) was multi- 
plied by 0.5 to facilitate comparison. The curve in (a) is a fit to the background distribution 
described in the text, while in (b) the shape from (a) is used with the normalization allowed 
to float and a signal Gaussian of width 2.7 MeV is added. 



We repeat this procedure for events in the oj sidebands. We use for our x\ definition pseudo-w 
masses in the sideband intervals. We show the Mb distribution for events in the AE sideband, 



denned above, and those having AE within 2a around zero in Fig. 17. We find no significant signal. 



7.3 B° D+U7T- Signal 

The same procedure followed for the D° final state is used for the D + final state. We show the 
candidate B mass distribution, Mb, for events in the AE side-band on Fig. |l~8|(a). The AE 
resolution is 18 MeV (a). This gives a good representation of the background in the signal region. 



We fit this distribution with a shape given in equation 15 



We next view the Mb distribution for events having AE within 2a around zero in Fig. |lq(b). 
This distribution is fit with a Gaussian signal function of width 2.7 MeV and the background 
function found above whose normalization is allowed to vary. We find 91±18 events in the signal 
peak. 

We repeat this procedure for events in the uj sidebands. We show the Mb distribution for both 
AE sidebands and AE within 2a around zero in Fig. 

There is no evidence of any signal in the u sideband plot, leading to the conclusion that the 
signal is associated purely with oj. 
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Figure 17: The B candidate mass spectra for the final state D°ujtt~ , with D° — > K~tt + and 
lo sidebands (a) for AE sidebands and (b) for AE consistent with zero. 

7.4 Branching Fractions 

We determine the branching ratios, shown in Table |||, by performing a Monte Carlo simulation of 
the efficiencies in the two modes. We use the current particle data group values for the relevant co, 
D+ and D° branching ratios of (88.8±0.7)% (w Tr+vr"^ ), (9.0±0.6)% (D + K--k + it + ) and 
(3.85±0.09)% (D° — ► K~tt + ) [HJ. The efficiencies listed in the table do not include these branching 
ratios g. 



Table 3: Branching Fractions for the Dam final state 



D° 


Decay Mode 


Fitted # of events 


Efficiency 


Branching Fraction (%) 




7T + 


88±14 


0.064 


0.41±0.07±0.04 


K~ 


7T + 7T + 


91±18 


0.046 


0.28±0.05±0.03 



The systematic error arises mainly from our lack of knowledge about the tracking and tt° 
efficiencies. We assign errors of ±2.2% on the efficiency of each charged track, and ±5.4% for the 
7r°. The error due to the background shape is evaluated in three ways. First of all, we change 
the background shape by varying the fitted parameters by la. This results in a change of ±5.0%. 
Secondly, we allow the shape, P2, to vary (the normalization, pi, was already allowed to vary). This 
results in 5.5% increase in the number of events. Finally, we choose a different background function 

back' (r) = p\r\J 1 — r 2 ^1 + + p%r 2 + p^r^j , (16) 

and repeat the fitting procedure. This results in a 1.0% decrease in the number of events. Taking 
a conservative estimate of the systematic error due to the background shape we arrive at ±5.5%. 
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Figure 18: The B candidate mass spectra for the final state D~ 
(a) for AE sidebands and (b) for AE consistent with zero. The vertical scale in (a) was 
multiplied by 0.5 to facilitate comparison. The curve in (a) is a fit to the background distri- 
bution described in the text, while in (b) the shape from (a) is used with the normalization 
allowed to float and a signal Gaussian of width 2.7 MeV is added. 



7.5 The utt System 

For all subsequent discussions we add the D° and D + final states together. We select sample of w's 
in the 7r + 7r _ 7r° mass window of 782±20 MeV using only combinations having r < 0.7 in the Dalitz 
plot. 

mass spectrum in the left-side plot. The solid histogram shows 
5.257 GeV) suitably normalized. The dotted 



In Fig. 20 we show the w-zr" 
events from the lower Mb sideband region (5.203 
histogram shows the background estimate from the AE sidebands, again normalized. In the signal 
distribution there is a wide structure around 1.4 GeV, that is inconsistent with background. We 
re-determine the (jtt~ mass distribution by fitting the Mb distribution in bins of lutt~ mass, and 
this is shown on the right-side. Fitting to a Breit-Wigner function, we find a peak value of 1415±43 
MeV and a width of 419±110 MeV. It should be kept in mind that this particular signal shape 
is assumed to be correct. Other resonant or non-resonant contributions could affect the mass and 
width. 

This structure appears identical to the one we observed in B — > D*ojit~ decays. 



7.6 Angular Distributions in Dun 

We can determine the spin and parity of the A~ particle by studying the angular distributions 
characterizing its decay products. The decay chain that we are considering is B — > A D; A — » ujtt 
and uj — > 7r + 7r~7r°. The helicity formalism || is generally used in the analysis of these sequential 
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Figure 19: The B candidate mass spectra for the final state D + ujtt~ , with D + 
and us sidebands (a) for AE sidebands and (b) for AE consistent with zero. 



K 7T + 7r + 



decays. This formalism is well suited to relativistic problems involving particles with spin S and 
momentum p because the helicity operator h = S ■ p is invariant under both rotations and boosts 
along p. 

There are two relevant reference frames. The first one, that we will define xaVaza is the rest 
frame of the A particle, with the z"a axis pointing in the A direction of motion in the B rest frame. 
The second one, x UJ y UJ z w , is related to xaVaza by the rotation through 3 Euler angles <pA, @Ai —4>A, 
as shown in Fig. ^lj. The angle (ft a defines the orientation of the plane containing the uj direction in 
the A rest frame and the za axis with respect to the x"a — z"a plane. The x"a direction is arbitrary. 
The angles W and (ft w define the orientation of the u decay plane in the uj rest frame. Note that 
the A decay plane has an azimuthal angle (ft>A both in the xaVaza and in the x w y w z w references. As 
the angle (ft a is arbitrary, the only angle that has a physical meaning is x = $A — (ftuj, the opening 
angle between the A decay plane and the to decay plane. 

Both the B meson and the D meson are pseudoscalar, therefore their helicity is 0. Thus A will 
be longitudinally polarized independently of its spin. In order to calculate the decay amplitude for 
this process, we need to sum over the uj helicity states: 



A 



(17) 



here D x l ((ft a, 9 a,— 4> a) is the rotation matrix that relates the xaVaza and the x LO y LO z w frames and 
D^ 1 {^un^u), —(ftu) is the rotation matrix relating the x^y^z^ and the direction of the normal to the 
uj decay plane n(6u,(ftu). 

In general, there are three helicity amplitudes that contribute to this decay: B\q and -B_io, cor- 
responding to a transverse uj polarization, and Bqq, corresponding to a longitudinal uj polarization. 
This expression can be simplified by observing that A — > ujtt is a strong decay and thus conserves 
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Figure 20: The invariant mass spectra of um~ for the final state Dujtc~ for both D decay 
modes, (left) The solid histogram is the background estimate from the Mb lower sideband 
and the dashed histogram is from the AE sidebands; both are normalized to the fitted 
number of background events, (right) The mass spectrum determined from fitting the Mb 
distribution and fit to a Breit-Wigner function. 



parity. Thus, the helicity amplitudes are related by the equation: 

B w = (-I^-^tmt^S-io, (18) 

#oo = {-lf- S(A) VA^B m . (19) 

where S(A) is the spin of particle A and an d represent the intrinsic parity of the decaying 

particle and its decay products, respectively. 



Eq. 18 relates the two transverse helicity amplitudes, while Eq. IS forbids the presence of a 
longitudinal component under certain conditions. For example, if A is a l - object, uj has transverse 
polarization and -B_io = —-Bio- When the sign in Eqs. 18-1^ is positive, two parameters determined 



by the hadronic matrix element affect the angular distribution and thus we cannot fully determine 
it only on the basis of our assumptions on the A spin parity. We have carried out the calculation of 
the predicted angular distributions including spin assignment for A up to 2. The predicted angular 
distributions are summarized in Table |||. 

The statistical accuracy of our data sample is not sufficient to do a simultaneous fit of the joint 
angular distributions shown above. Thus only the projections along the 9a, @u and \ are fitted, 
integrating out the remaining degrees of freedom. Table [|, gives the analytical form for these 
projections. 

We determine the projections of these angular distributions by fitting the Mb distribution as 
a function of the various angular quantities cosOa, cos# w , %■ We restrict the u>ir~ mass range 
to be between 1.1 and 1.7 GeV, containing 104 events. In order to fit the angular distribution 
with theoretical expectations, we must the correct the data for acceptances. We determine the 
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Figure 21: Relationship between the A rest frame xaUa z a and the u rest frame x^y^z^. xa 
and x w lie in the same plane. 



Table 4: Differential angular distributions (modulo a proportionality constant) predicted for 
different spin assignments. 



J v da / d cos 6 Ad cos 8 ^dx 

0^ \B 00 \ 2 cos 2 9 U 

1 ~ | B w | 2 sin 2 9 a sin 2 9 U sin 2 x 

1 + | Bio 1 2 sin 2 6U sin 2 B u 2 cos x 2 + | Boo | 2 cos 2 9 A cos 2 W 

-l/2Be(B 10 B* ) sin 20 A sin 20 w cos x 
2" 3\B W \ 2 sin 2 20 A sin 2 9 W cos 2 x + |B 00 | 2 (3 cos2 9 a ~ I) 2 cos 2 ( 

-V3Re(B 10 B* ) sin 20 A (3 cos 2 A - 1) sin 29 w cos % 
2+ 3/4 1 B 10 1 2 sin 2 20 A sin 2 W sin 2 x 



Table 5: Projection of the angular distributions along the cos 0,4, cos9 UJ and x axes. 



J v da /d cos 9 a da /d cos 9^ da/dx 

0^ f\B^ 2 47r|Boo| 2 cos 2 w 4/3|B 00 | 2 

1- f |B lo | 2 sin 2 yl f |B lo | 2 sin 2 w 1 1 B w | 2 sin 2 X 
1+ f (|B lo | 2 sin 2 A f(|B lo | 2 sin 2 w §(4|B 10 | 2 cos 2 X 

+\B 00 \ 2 cos 2 9 A ) + |B oo | 2 cos 2 0J + |B 00 | 2 ) 

2- f (3|B lo | 2 sin 2 20 A ^(|B 10 | 2 sin 2 W 4|B 10 | 2 cos 2 X 
+ |B oo | 2 (3cos 2 A -l) 2 ) + |B oo | 2 cos 2 w ) + |B 00 | 2 

2+ 7r|B lo | 2 sin 2 20 A f |B 10 | 2 sin 2 0^ jf|B 10 | 2 sin 2 x 
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acceptance correction by comparing the Monte Carlo generated angular distributions with the 
reconstructed distributions. The angular dependent efficiencies are shown in Fig. 22. 
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Figure 22: Reconstruction efficiency dependence on (a) cos 9a, (b) cos 6^, and (c) x- 



The corrected angular distributions are shown in Fig. 23. The data are fit to the expectations 
for the various J p assignments. For the CP, 1~ and 2 + assignments, the curves have a fixed shape. 
For the 1 + and 2~ assignments we let the ratio between the longitudinal and transverse amplitudes 
vary to best fit the data. We notice that the u> polarization is very clearly transverse (sin 2 9 W ) and 
that infers a 1~ or 2 + assignment. 

We list in Table |6| the x 2 /dof for the different J p assignments. The 1~ assignment is preferred, 
having a x 2 /dof of 1.7. The other assignments are clearly ruled out. The probability that we have 
a correct solution and x 2 /dof is 1.7 or greater is 3.8% @. 



Table 6: The x 2 of angular fits 





cr 


1+ 


1- 2+ 


2- 


X 2 /dof 


7.0 


4.5 


1.7 3.2 


5.3 


probability 


1.9 x 1CT 15 


3.3 x 1CT 8 


3.8% 2.7 x 1CT 5 


3.3 x TCP 10 



8 Discussion of Nature of the A 

We have found a l - object decaying into ujtt~ . A simple Breit-Wigner fit assuming a single 
resonance and no background gives a mass of 1418±26 MeV with an intrinsic width of 382±41 
MeV. (The individual measurements are listed in Table 0). We have evaluated the systematic 
errors due to changes in the parameterization of the background shape in when fitting the Mb 
distributions in bins of luit~ mass. This leads to an error in the mass of 19 MeV and 32 MeV in 
the width. 
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Figure 23: The angular distribution of 9 A (top-left), 6 U (top-right) and x (bottom). The 
curves show the best fits to the data for for different J p assignments. (The 0~ and 1 + are 
almost indistinguishable in cos #4, while the 1~ and 2 + are indistinguishable in cos6^ and 

X- 
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Table 7: Measured p' Mass and Width 



Mode 


Mass (MeV) 


Width (MeV) 




1432±37 


376±47 




1367±75 


439±135 




1415±43 


419±110 


Average 


1418±26 


388±41 



Other possible changes in the mass and width are more difficult to evaluate. They include 
different Breit-Wigner parameterizations and the possibility of other resonant or non-resonant 
components in the ujtt~ mass distributions. 

Signals for um~ resonances have been detected before below 1500 MeV. There is a well es- 
tablished axial- vector state the bi(1235) with mass 1230 MeV and width 142 MeV. Data on vec- 
tor states, excited p's, is inconsistent. Clegg and Donnachie [11| have reviewed — ► (4-7r)~P, 
e + e~ — > ir + ir~ and e + e~ — > 7r + 7r + 7r~7r~ data, including the uitt final state. Their best explanation 
is that of two 1" states at 1463±25 MeV and 1730±30 MeV with widths 311±62 and 400±100 
MeV, respectively. Only the lighter one decays into coir. The situation is quite complex, however. 
They conclude that these states must be mixed with non-gg states in order to explain their decays 
widths. There is also an observation of a wide, 300 MeV, ujtt° state in photoproduction at 1250 
MeV g|], that is dominantly the &i(1235) @ with possibly some 1 in addition. Our state is 
consistent with the lower mass p' . We do not seem to be seeing significant production of the higher 
mass state into luit~ , as expected. 

Several models predict the mass and decay widths of excited p and u> mesons . For example, 
according to Godfrey and Isgur [14| the first radial excitation of the p is at 1450 MeV. There is 
a large variation among the models, however, on prediction of the relative decays widths ranging 
from no tttt to tttt being equal to loit [15]. 

Since we have observed a wide 1~ state in the mass region where the p' is expected, the most 
natural explanation is that we are observing the p' for the first time in B decays. 

We note that t~ lepton decays into ujit~ have been observed, and the 1~ spin-parity defi- 
nitely established [|l(|. However, the relatively low mass of the r~ distorts the the mass spectrum 
significantly, and makes it difficult to extract the p' mass and width [17]. 



9 Search For Other Resonant Substructure in D*(4ir) 

We have accounted for ~20% of the (47r)~ final state. We would like to disentangle other resonant 
substructure. Since the background is large in modes other than D° — > K~ir + we will only use this 
mode. One process that comes to mind is that where the virtual W~ materializes as an oj~, that 
subsequently decays into ir + ir~ir~ and we produce a D** + that decays into D* + ir°. This process 
should be the similar to that previously seen in the reaction B~ — * D**°ir~ , where the D**° decayed 
into a D* + n~ [jlq l. We search for the presence of an by examining the 7r + 7r _ 7r _ mass spectrum 



in Fig. 24 



There is an excess of signal events above background in the aj - mass region, that cannot be 
definitely associated with the a±. Proceeding by selecting events with ir + ir~ir~ masses between 0.6 
and 1.6 GeV, we show the D* + ir° invariant mass spectrum in Fig. p^. 
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Figure 24: The invariant mass spectra of 7r + 7r~7r~ for the final state D* + n + n 
D° — > K~7i + . The solid histogram is the background estimate from the Mb lower sideband 
and the dashed histogram is from the AE sidebands; both are normalized to the fitted 
number of background events. 
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Figure 25: The invariant mass spectra of D* + tt° for 7r + 7r~7r~ masses between 0.6 - 1.6 GeV 
for the final state D* + n + 7Y^n^n with D° — > K~it + . The solid histogram is the background 
estimate from the M B lower sideband and the dashed histogram is from the AE sidebands; 
both are normalized to the fitted number of background events. 
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Although there is a suggestion of a low mass enhancement, it is not consistent with D** pro- 
duction that would peak in region of 2.42 - 2.46 GeV. Perhaps we are seeing an indication of 
fragmentation at the b — > c decay vertex here. 

We also display for completeness the "aj~7r°" mass distribution in Fig. 26. There may or may 
not be a wide structure in the (47r)~ mass. At this point we abandon our search for substructure 
in this decay channel. 
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Figure 26: The invariant mass spectra of 7r + 7r~7r~7r° for 7r + 7r~7r - masses between 0.6 - 
1.6 GeV for the final state D* + 7t + 7t~7t~tt with D° — ► K~n + . The solid histogram is the 
background estimate from the Mb lower sideband and the dashed histogram is from the AE 
sidebands; both are normalized to the fitted number of background events. 



10 Conclusions 

We have made the first statistically significant observations of six hadronic B decays shown in 
Table |. 

There is a low-mass resonant substructure in the ujit~ mass. A simple Breit-Wigner fit assuming 
a single resonance and no background gives a mass of 1418±26±19 MeV with an intrinsic width of 
382±41±32 MeV. 

The structure at 1418 MeV has a spin-parity consistent with 1~. It is likely to be the elusive 
p' resonance [PH| . These are by far the most accurate and least model dependent measurements of 
the p 1 parameters. The p' dominates the final state. (Thus the branching ratios for the D^ujtt^ 
apply also for D^p'~.) 

Heavy quark symmetry predicts equal partial widths for D* p' and Dp'. We measure the relative 
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Table 8: Measured Branching Ratios 



Mode 




B (%) 


of events 


B u -> 


D* + 7l + 7l-7l-7T° 


1.72±0.14±0.24 


1230±70 


B° 


D* + uotj- 


0.29±0.03±0.04 


136±15 


B° 


D + U7T- 


0.28±0.05±0.03 


91±18 


B~ - 


D*°TT + TT-TT-TT° 


1.80±0.24±0.25 


195±26 


B~ - 


D*°LUTT- 


0.45±0.10±0.07 


26±6 


B~ - 


D°luti- 


0.41±0.07±0.04 


88±14 



rates to be 

r (b° -> d*+ p '- 



r (b° -> L>V 

r (5- -> 
r(#- oy 1 ; 
r(s -> d* p '~) 



■ = 1.04 ±0.21 ±0.06 (20) 

= 1.10 ±0.31 ±0.06 (21) 
1.06 ±0.17 ±0.04 . (22) 



T(B^Dp'-] 

Thus the prediction of heavy quark symmetry is satisfied within our errors. 

Factorization predicts that the fraction of longitudinal polarization of the D* + is the same as in 
the related semileptonic decay B — > D*l~v at four-momentum transfer q 2 equal to the mass-squared 
of the p' 

Y L {B -» D*+p'-) _T L (B^ D*£~v) 



T(B^D*+p'-) T(B^D* 



(23) 



Our measurement of the D* + polarization (see Fig. 14) is (63±9)%. The model predictions in 
semileptonic decays for a q 2 of 2 GeV 2 , are between 66.9 and 72.6% |jl9f . Thus this prediction of 
factorization is satisfied. 

We can use factorization to estimate the product of the p' decay constant f p i and the branching 
ratio for p'~ — > lott~ . The relevant expression is 



£ (b - nn-v) \ q2 =m * 



6TT 2 clflB{p'~ ^LUir-)\V ud \ 2 , (24) 



where f p > is the so called p' decay constant and c\ is a QCD correction factor. We use ci=l.l±0.1 

We use the semileptonic decay rates given in Barish et al. [21]. The product 

f 2 ,B (p'~ um~) = 0.011 ± 0.003 GeV 2 , (25) 

where the error is the quadrature of the experimental errors on the experimental branching ratios 
and c\. 

The model of Godfrey and Isgur predicts decay constants widths and partial widths of mesons 
comprised of light quarks by using a relativistic treatment in the context of QCD 0]. They predict 
both f p i and B (p'~ — > ujtt~); the values are 80 MeV and 39%, respectively. The branching ratio 
prediction is believed to be more accurate (p2|. We use this to extract 



fpf = 167 ± 23 MeV . (26) 
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The model predicts a lower value for f p i observed here, if factorization is correct. 

We note that all the B — > D^p' branching ratios that we have measured are approximately 
equal to the B — > D^p branching rates [Q] if a model value of B (p'~ —> ujit^) = 39% is used. 

Finally, although the B° — > D* + (4ir)~ and B~ —> D*°(<iir)~ branching ratios are nearly equal, 
the u>ir~ branching ratios are about 1.5 times larger for the charged B than the neutral B, main- 
taining the trend seen for the ir~ and p~ final states. Since the B~~ lifetime is if anything longer 
than the B°, this trend must reverse for some final states. It has not for D^p'. 
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